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ABSTRACT 

In Vanderbei et al. (2003), we introduced a class of shaped pupil masks, called 
spiderweb masks, that produce point spread functions having annular dark zones. 
With such masks, a single image can be used to probe a star for extrasolar plan- 
ets. In this paper, we introduce a new class of shaped pupil masks that also 
provide annular dark zones. We call these masks starshaped masks. Given any 
circularly symmetric apodization function, we show how to construct a corre- 
sponding starshaped mask that has the same point-spread function (out to any 
given outer working distance) as obtained by the apodization. 

Subject headings: Extrasolar planets, coronagraphy, Fraunhaufer optics, Hankel 
transform, point spread function, apodization 



1. Introduction 



With more than 100 extrasolar Jupiter-sized planets discovered in just the last decade, 
there is now great interest in discovering and characterizing Earthlike planets. To this end, 
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NASA is planning to launch a space-based telescope, called the Terrestrial Planet Finder 
(TPF), sometime in the middle of the next decade. This telescope, which will ultimately 
be either an interferometer or a coronagraph, will be specifically designed for high-contrast 
imaging. Earlier studies (Brown et al. (2002)) indicate that a 4m class coronagraph ought 
to be able to discover about 50 extrasolar Earth-like planets if the telescope can provide 
contrast of 10~ 10 at a separation of 3X/D and that a 4 x 10m class telescope ought to be 
able to discover about 150 such planets if it can provide the same contrast at a separation 
of AX/D. 

One of the most promising design concepts for high-contrast imaging is to use pupil 
masks for diffraction control. Early work in this direction (see Brown et al. (2002); Spergel 
(2000); Kasdin et al. (2002, 2003)) has focused on optimizing masks that are not rotationally 
symmetric and thus provide the desired contrast only in a narrow annular sector around a 
star, but at fairly high throughput. A full investigation of a given star then requires multiple 
images where the mask is rotated between the images so as to image all around the star. In 
a recent paper (Vanderbei et al. (2003)), we proposed a new class of rotationally symmetric 
pupil masks, which are attractive because they do not require rotation to image around 
a star. Such masks consist of concentric rings supported by many ultra-thin pie-shaped 
spiders. We call these masks spiderweb masks. We showed that they can provide the desired 
contrast with a reasonable amount of light throughput. 

In this paper, we consider a similar family of binary masks. These masks consist only of 
spiders — there are no concentric rings. Instead, the desired contrast is obtained by carefully 
controlling the width of the spiders as a function of radius. We call these masks starshaped 
masks since the resulting designs look like many-pointed stars. 

The paper is organized as follows. In the next section, we briefly review the relationship 
between pupil-plane apodization and the corresponding image-plane electric field and point- 
spread function. Then in Section 3, we show how to make starshape binary masks for 
which the first term in a Bessel expansion of the electric field matches the corresponding 
electric field associated with any given apodization. We show that for a sufficiently large 
number of star-points, the remaining terms have negligible contribution to the electric field 
in the annulus of interest. Finally, in Section 4 we show how to compute various "optimal" 
apodizations. In particular, we show that the best of these optimal apodizations is zero-one 
valued and in fact corresponds to the optimal spiderweb masks presented in Vanderbei et al. 
(2003). 
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2. Apodization 



The image-plane electric field produced by an on-axis point source and an apodized 
aperture defined by a circularly-symmetric apodization function A(yJ x 2 + y 2 ) is given by 



and r(x,y) and 9(x,y) denote the polar coordinates associated with point (x,y). Here, and 
throughout the paper, x and y denote coordinates in the pupil plane measured in units of the 
aperature D and £ and ( denote angular (radian) deviations from on-axis measured in units 
of wavelength over aperture (X/ D) or, equivalently, physical distance in the image plane 
measured in units of focal-length times wavelength over aperture (fX/D). 

For circularly-symmetric apodizations, it is convenient to work in polar coordinates. To 
this end, let r and 9 denote polar coordinates in the pupil plane and let p and denote the 
image plane coordinates: 





s 



where 



S={(x,y): < r(x, y) < 1/2, 9(x, y) G [0, 2vr]} , 



(2) 



x = r cos 9 
y = r sin 9 



£ = p cos 
C = psin^. 



(3) 



Hence, 



rp(cos 9 cos + sin 9 sin 0) 
rpcos(9 — 0). 



(4) 
(5) 



The electric field in polar coordinates depends only on p and is given by 



1/2 2tt 




(6) 







1/2 




(7) 







where Jo denotes the 0-th order Bessel function of the first kind. Note that the mapping 
from apodization function A to electric field E is linear. Furthermore, the electric field in 
the image plane is real- valued (because of symmetry). 



-4- 



The point spread function (psf) is the square of the electric field in the image plane. The 
contrast requirement is that the psf in the dark region be 10~ 10 of what it is at its center. 
Because the electric field is real-valued, it is convenient to express the contrast requirement 
in terms of it rather than the psf, resulting in a field requirement of ±10~ 5 . 

A second requirement is high throughput. The most natural measure of throughput is 
the area under the main lobe of the psf: 

Piwd 

T Airy = J E 2 (p)2irpdp, (8) 
o 

where p iw d denotes the location of the first null of the psf. We call this measure the Airy- 
throughput. Two other measures of throughput are relevant to our discussion. The total 
throughput of an apodization is the integral of the psf over the entire image plane. By 
Parseval's theorem, the total throughput is also the integral of the square of the apodization: 

oo 1/2 

T tot&l = J E 2 {p)2npdp = J A 2 {r)2itrdr. (9) 
o o 

The electric field at p — 0, 

1/2 

E(0) = J A(r)27rrdr, (10) 
o 

provides another measure of the "central" throughput of the apodization, since its square is 
the peak throughput density at the center of the Airy disk. If AQ is zero-one valued, then 
the apodization can be realized as a mask. In this case, E(0) is precisely the open area of 
the mask (and is also the total-throughput). For this reason, we call E(0) the pseudo-area 
of the apodization even when the apodization is not zero-one valued. 



3. Starshape Masks 



In this section, we study binary mask approximations to any apodized pupil. The masks 
we consider are star shaped — see Figure 1. The opening in an Appoint starshaped mask can 
be described mathematically by a set S given by: 



S = {(x,y): < r(x, y) < 1/2, 0(x, y) G 0} , 

N-l 



n=0 



2nn a(r) 2ir(n + 1) a(r) 



N 



N 



(11) 
(12) 
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where a(r) denotes the width in radians of a "vane" and the notation [a, b] denotes the 
interval on the real line from a to b. Note that the shape of each point of the star is 
determined by the function a(r). Our aim is to determine choices for this function that will 
yield an image-plane psf matching the psf corresponding to any given apodization. 

The electric field, expressed in polar coordinates, associated with starshape mask S is 
given by 

E{p,(j>) = J J e- 2mrpcos{e -^rdrd6. (13) 
s 

The integral in equation (13) can be expressed in terms of Bessel functions using the 
Jacobi- Anger expansion (see, e.g., Arfken and Weber (2000) p. 681): 

oo 

e lxcose = i m J m {x)e ime . (14) 

m=— oo 

Substituting into (13), we get: 

E{p,</>) = I j^2i m J m {-2Tirp)e im ^ ) rdrde (15) 

s m 

= ! ^i m J m {-2Tirp)e- im ^ j f e ime de \ rdr (16) 
o m \e / 

The integral over 6 is easy to compute: 

2tt(ti+1) a(r) 
JV-1 N 2~ 

j e ime d9 = J eim6d9 ( 17 ) 



e n ~° 2nn , a(r) 

N 2 



2tt - Na(r) m = 

-fsin(jA^r)) m = jNij ^ Q (18) 
otherwise. 



Substituting this result into (13), yields 

1/2 



E(p,(f) = J J (-2irrp)(2n - Na{r))rdr (19) 
o 

1/2 

- / i jN JjN(- 2 ^p) e ~ ijN(t> - sin (j N ^r) rdr - ( 2 °) 
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Lastly, suppose that iV is even and use the fact that J- m (x) = J m (—x) = (—l) m J m (x) to 
get the following expansion for the electric field: 



1/2 

f N 

E(p,(j)) = 2ir / J (27rrp)(l-— a(r))rdr 



(21) 



oo 1/2 

4 5^ / ^(27rrp)cosO'7V(0-7r/2))-sinOW^i)rdr. (22) 



The first term, involving the integral of Jo, is identical to the formula for the electric 
field for an apodized aperture with 



Hence, by putting 



N 

Mr) = 1 - ^a(r) 



2tt 

a(r) = -(l-A(r)), 



(23) 



(24) 



we can make the first term match any circularly symmetric apodization. 

Furthermore, for large N, the effect of the higher-order Bessel terms becomes negligible 
for small p. Indeed, for z < ^4(m + 1), 



< J m (z) < 



(*/2) 



m+1 



(m + 1)! ' 



(25) 



(which itself follows easily from the alternating Taylor series expansion of the m-th Bessel 
function: J m (z) = X!z^o( — ^-) l (z /2) 2l+m / {l\{m+l)\)). From this we use the Schwarz inequality 
to estimate the magnitude of the effect of the higher-order terms: 



/ J jN (2nrp)cos(jN(<f)-n/2))-sm(jN^p-)rdr 
7-i •* J 2 

J— 1 o 

1 x 

-2^(jN + l)\\ 2 J 



(26) 



for p ow( j < a/4(A^ + l)/vr. Since the last bound is dominated by e^owd/ 2 /2, it follows from 
the dominated convergence theorem that this last bound tends to zero as N tends to infinity. 
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The convergence to zero of the terms in the sum on j is very fast. In fact, if N is set 
large enough such that max < z < 7rPowd Jn(z) < 10~ 5 , then the j = 1 term dominates the sum 
of all the higher-order terms and is itself dominated by the Jo term. Figure 2 shows a plot of 
J 50 , Jioo> an d Ji5o- These three Bessel functions first reach 10" 5 at z = 35.2, 81.0, and 128.1, 
respectively. This suggests that a contrast level of 10~ 10 can be preserved out to 35X/D 
using about N = 50 vanes, out to 80X/D using about N = 100 vanes, and out to 120X/D 
using about N = 150 vanes. These estimates ignore the contribution from higher-order 
terms. In practice, more vanes are required to preserve the desired contrast level to these 
outer working distances. 

Within the discovery zone, the electric fields resulting from the apodization and the 
starshaped mask agree. Hence, the Airy-thoughputs and pseudo-areas also agree. But, the 
total throughput for the mask is larger, often by a factor of about two. Indeed, the total 
throughput of the apodization is 

1/2 

r tota i, apod = J A(r) 2 2nrdr (27) 
o 

whereas, from (12) and (23), it follows that the total throughput of the mask is 

1/2 

^otai.mask = J J rdrdO = J A(r)2nrdr > T apod (28) 

S 

(the inequality holds because < A(r) < 1). So, the mask lets about twice as much light 
through but the extra light is necessarily concentrated outside the outer working distance. 



4. Apodizations Optimized for TPF 

There has been much interest in using apodized pupils to achieve high-contrast imaging, 
especially recently in the context of TPF studies: see e.g. Jacquinot and Roizen-Dossier 
(1964); Indebetouw (1990); Watson et al. (1991); Nisenson and Papaliolios (2001). 

The simplest optimization problem involves maximizing the pseudo-area subject to con- 
trast constraints. It can be formulated as an infinite dimensional linear optimization problem: 

1 /9 

maximize 

subject to -10~ 5 £(0) < E(p) < 10- 5 £(0), p iwd < p < p owd , (29) 

0<A(r)<l, 0<r<l/2, 
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where p; w d denotes a fixed inner working distance and p owc j a fixed outer working distance. 
Discretizing the sets of r's and p's and replacing the integrals with their Riemann sums, 
problem (29) is approximated by a finite dimensional linear programming problem, which 
can be solved to a high level of precision (see, e.g., Vanderbei (2001)). The numerical solution 
to this problem reveals that the optimal solution is zero-one valued. That is, the optimal 
apodization is, in fact, a mask consisting of concentric rings. 

A second optimization problem involves maximizing total-throughput: 



An empirical observation is that this more difficult optimization problem also proves to be 
tractible and, in all cases tried, provides the same optimal apodization as was obtained by 
maximizing pseudo-area. 

The solution obtained for p iw< \ = 4 and p ow d = 60 is shown in Figure 3. The throughput 
is 17.9% and, since it is a mask, its pseudo-area is the same. Its Airy-thoughput is 9.37%. 
(Note: all results are reported as a percentage of the largest possible value which would 
result from a completely open aperture. That is, all throughputs are divided by the area 
7r(l/2) 2 and labeled as percentages.) 

The starshaped mask associated with this "apodization" is precisely the same concentric- 
ring mask independent of the number of star points chosen, since the star points disappear 
entirely. Of course, this mask cannot be manufactured as the concentric rings need some 
means of support. If one imposes an upper bound on A in (29) that is strictly less than 
unity, then one obtains a supportable mask with slightly reduced throughput (the amount 
depending on the upper bound chosen). Such a mask is exactly the spiderweb mask studied 
in Vanderbei et al. (2003). 

In this paper, we are interested in smooth apodizations. To this end, we add seemingly 
artificial smoothness constraints to the optimization problem. Motivated by the fact that 
optimal apodizations look qualitatively like a Gaussian function, we impose smoothness 
constraints that correspond to the following conditions: 



maximize J L/2 A(r) 2 27rrdr 
subject to -10- 5 £(0) < E(p) < 10- 5 £(0), 

< Mr) < 1, 



Piwd < P < Powd, 

< r < 1/2. 



(30) 



log(A)' < 
log(A)" < 0. 



(31) 
(32) 
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The resulting max-pseudo-area optimization problem is: 



maximize J Q A(r)27rrdr 

subject to -l(r 5 £(0) < E(p) < 10- 5 £(0), 

< A(r) < 1, 
A'{r) < 0, 
A(r)A"(r) < A'(r) 2 , 



Piwd < P < Powd, 



< r < 1/2, 
< r < 1/2, 
< r < 1/2. 



(33) 



The solution obtained for p iw< j = 4 and p ow d = 60 is shown in Figure 4. The total- 
throughput for this apodization is 9.12%. Its Airy-throughput is 9.09% and its pseudo-area 
is 17.39%. The psf's for the corresponding starshape masks with 20 and 150 point stars 
are shown in Figure 5. Other apodizations satisfying the smoothness constraints (such as 
the generalized prolate spheroidal apodization discussed in the next subsection), being less 
optimal, must necessarily have even smaller throughputs. 

There is great interest in reducing the inner working distance to something less than 
Piwd = 4, as this would increase the number of likely target stars for TPF. Smaller inner 
working distances can be achieved if if one is willing to accept a smaller outer working distance 
too (di Francia (1952)). Unfortunately, a small reduction in inner working distance requires 
a large reduction in outer working distance so that very quickly the discovery zone becomes 
a very narrow annulus. For example, the solution obtained for pi w< j = 3 and p owc j = 4.25 is 
shown in Figure 6. The total-throughput for this apodization is 7.7%. Its Airy-throughput 
is 7.6% and its pseudo-area is 20.0%. The 50-point starshaped mask associated with this 
apodization is shown in Figure 7. 

Optimizing pseudo-area or total-throughput are really just simple surrogates for the 
true objective, which is to minimize integration time. Burrows (1994) has shown that maxi- 
mizing sharpness translates directly into minimizing the integration time required to reach a 
specified signal to noise ratio. However, it seems to us that optimizing sharpness is beyond 
the capabilities of our optimization tools. Hence, we optimize the simple surrogates and 
reserve sharpness calculations for later comparing various designs. 



Generalized prolate spheroidal wave functions were first introduced by Slepian (1965) 
as a way to apodize a circularly symmetric aperture so as to concentrate as much light as 
possible into a central Airy disk. In the context of TPF, the generalized prolate spheroidal 
wave function has been popularly viewed as providing an optimal apodization (see, e.g., 
Kasdin et al. (2003); Aime et al. (2001); Gonsalves and Nisenson (2003)). However, in truth, 



4.1. The Generalized Prolate Spheroidal Apodization 
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it is an optimal solution to the following slightly different optimization problem 1 : 



minimize 




(34) 



Piwd 



subject to 



A(0) = 1. 



The generalized prolate spheroidal apodization computed using p iwd = 4 gives 8.3% total- 
throughput. The Airy-throughput is essentially the same. The pseudo-area is 15.9%. The 
apodization and corresponding psf are shown in Figure 8. At the expense of Airy-throughput 
(8.3% vs. 9.1%) and slight violation of the contrast in the first few diffraction rings, it has 
better than needed contrast throughout most of a larger than needed dark zone. 

The throughputs for the apodizations discussed in this section are summarized in Table 

1. 



Some TPF concepts involve an elliptical pupil geometry since this might provide a means 
to achieve improved angular resolution in realizable rocket fairings. The designs presented 
in this paper are given in unitless variables. When re-unitizing, a different scale can be used 
for the x and y directions. In this way, these designs can be applied directly to elliptical 
pupils. Of course, the high-contrast region of the psf will also be elliptical with the short 
axis of the psf corresponding to the long axis of the pupil. 

There are two parts to TPF: discovery and characterization. Discovery refers to the 
simple act of looking for exosolar planets. Characterization refers to the process of learning 
as much as possible about specific planets after they have been discovered. The masks 
presented here are intended primarily for discovery since a single exposure with these masks 
can discover a planet in any orientation relative to its star. However, once a planet is 
found and its orientation is known, some of the asymmetric masks presented in previous 
papers should be used for photometry and spectroscopy as they have significantly higher 
single-exposure throughput. 

An alternative to pupil masks is to use a traditional coronagraph, which consists of an 
image plane mask followed by a Lyot stop in a reimaged pupil plane. Recently, Kuchner 



Slepian actually formulated it as an equivalent unconstrained optimization problem: 



5. 



Final Remarks 



mzxj^EipYdp/ j™E{p)Hp. 
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and Traub (2002) have developed band-limited image-plane masks that achieve the desired 
contrast to within 3X/D. A potential drawback of this approach is its sensitivity to pointing 
accuracy. Nonetheless, it appears very promising. In the future, we plan to consider combin- 
ing pupil masks with image masks and Lyot stops to make a hybrid design that hopefully will 
provide a design achieving the desired contrast with an even smaller inner working distance 
and perhaps higher throughput. 

In this paper we have only considered scalar electric fields. We leave the important and 
more complex issue of how to treat vector electric fields, i.e. polarized light, to future work. 
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Fig. 1. — A 20-point starshape mask corresponding to the apodization shown in Figure 4. 




Fig. 2. — The Bessel functions J 50 , Jioo, and Jiso- They first reach 10 5 at 35.2, 81.0, and 
128.1, respectively. 
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Fig. 3. — Top. A concentric-ring mask designed to provide high-contrast, 10~ 10 , from X/D = 
4 to X/D = 60. Total throughput and pseudo-area are 17.9%. Airy throughput is 9.37%. 
Bottom. The associated psf. 
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Fig. 4. — The optimal apodization for p iw a = 4 and p ow d = 60 and the associated psf. Total 
throughput is 9.12%. Pseudo-area is 17.39%. Airy-throughput is 9.09%. 




Fig. 5. — Psf's for starshape mask in associated with apodization shown in Figure 4. Top 
Row. 20-point star. Second Row. 150-point star. 
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Fig. 6. — The optimal apodization for p iwc j = 3 and p ow< j = 4.25 and the associated psf. 
Total throughput is 7.7%. Pseudo-area is 20.0%. Airy-throughput is 7.6%. 




Fig. 7. — Psf 's for starshape mask in associated with apodization shown in Figure 6. 50-point 
star. 
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Fig. 8. — Left. The generalized prolate spheroidal apodization computed using p iw a = 4. 
Total throughput and Airy-throughput agree and are equal to 8.52%. The pseudo-area is 
16.28%. Right. The associated psf. 
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Smooth ? 


Piwd 


Powd 


^total 


^Airy 


E(0) 


Figure 3 


No 


4 


60 


17.9 


9.37 


17.9 


Figure 4 


Yes 


4 


60 


9.12 


9.09 


17.39 


Figure 6 


Yes 


3 


4.25 


7.7 


7.6 


20.0 


Figure 8 


Yes 


4 


60 


8.34 


8.34 


15.94 



Table 1: Summary of the throughput results for the specific apodizations considered. 



